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Abstract 

We study the conditions under which the nucleons inside a de- 
formed nucleus can undergo chaotic motion. To do this we per- 
form self-consistent calculations in semiclassical approximation utiliz- 
ing a multipole-multipole interaction of the Bohr-Mottelson type for 
quadrupole and octupole deformations. For the case of harmonic and 
non-harmonic static potentials, we find that both multipole deforma- 
tions lead to regular motion of the collective coordinate, the multipole 
moment of deformation. However, despite this regular collective mo- 
tion, we observe chaotic single particle dynamics. 



PACS: 21.10.Ky, 21.10.Re, 21.60.-n, 24.60.Lz 



The question about the origin of dissipation in collective motion of fi- 
nite Fermi systems [0 such as atomic nuclei or small metallic clusters is an 
intriguing and up to now not completely satisfactorily solved problem. For 
example, the mutual balance of one-body and two-body processes is still 
a question of debate. For the case of one-body dissipation and friction in 
nuclear dynamics, Swiatecki and coworkers 0, 0, |^ have developed this 
picture: particles which move in a shape-deformed container are reflected 
from the (moving) walls, and due to parts of it having positive curvature (for 
higher multipole moments) the particles very quickly loose their coherence, 
thus inducing pseudo-random motion, i.e. heat, into the system. At the same 
time the shape oscillation is very much slowed down. 

Blocki et al. 0, |[ consider a purely classical gas of particles contained 
in a deformed billiard. The only similarity with a Fermi gas comes from 
the fact that initially the particles' momenta are distributed within a Fermi 
sphere. The walls of the container undergo periodic shape oscillations with 
a frequency much smaller than a typical single particle frequency. In the 
interior of the container the particles move on linear trajectories. They study 
the particle kinetic energy increase as a function of time and find that for 
ellipsoidal shape deformations (£ = 2) the particles act as a classical Knudsen 
gas i.e. the total kinetic energy increase over an entire shape oscillation 
period is 0. However, for £ > 3 the kinetic energy in the single particle 
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motion is not completely 'given back', but rather steadily increases in time. 
This is explained by the fact that in an £ = 2 potential the motion of the 
particles remains non-chaotic and therefore coherent, whereas in the i > 3 
the scattering of the segments of the wall with positive curvature leads to 
chaotic motion similar to the one observed in a Sinai J^, H billiard and thus 
a destruction of coherence. The fact that deformed nuclear potentials may 
exhibit chaotic motion was recognized early by Arvieu and co-workers ^ 

This scenario is very similar to the so-called Fermi acceleration, proposed 



to explain the occurrence of very high energy cosmic radiation |10, 11 



In this paper we present an attempt to include selfconsistency into the 
problem of motion in multipole-deformed nuclear potentials. We have cho- 
sen a self consistent, but schematic, model of separable forces. We chose 
an interaction of the Bohr-Mottelson type |T2| with a static potential 



and multipole-multipole interactions as studied, for example, by Stringari et 
al. |l^, |l^, |15|, In the small amplitude limit, this model has recently 



been investigated in the semiclassical limit [|T6|; as is known, the low-lying 
quadrupole and octupole frequencies come out to be in reasonable agreement 
with experimental data. Our single-particle Hamiltonian is then 

= ^ + Vo + V^'\r,t) (1) 
2 m 
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V^^\r,t) is the potential associated with the (separable) multipole-multipole 
force [|13|, |I5 

V^'\r,t) = Xeqeir)Q,it) , (2) 

and Vq is the static external potential. We take here Vq = ^ m Uq r'^ , resulting 
in the Bohr-Mottelson Hamiltonian [|l^ . However, we have also investigated 
non-harmonic static potentials (r^) and obtained similar results. 

The coupling constants can be calculated using a self-consistent nor- 



malization condition |jT2|, [16 1. 



where A is the mass number of the nucleus under consideration. qe{r) is 
given by 

q2ir) = TyT^ (5) 

Q'iir) = r^TyT^ (6) 
and the multipole moments Quit) are 

f{r,p,t) is the one-body phase space distribution function of nucleons, the 
Wigner transform of the one-body density. 



We treat this problem in semi-classical approximation by a Wigner trans- 
formation of the von Neumann equation of motion for the density matrix, 
idtp = [H,p], to obtain a Vlasov equation, dtf = {H,f}. We then solve 
the Vlasov equation in the test particle method [l^, |l8| using a fourth-order 
Runge-Kutta algorithm with typical time step sizes of 1 fm/c. Our numerical 
simulation is fully selfconsistent and conserves total energy to better than 0.1 
%. 

In order to generate selfconsistent initial deformations in coordinate and 
momentum space, we start with a spherically symmetric configuration gener- 
ated in local Thomas- Fermi approximation without the deformation potential 

We now apply a time-dependent external potential of the form 

= kq sin(a'_D t) s{t) qe{r) (8) 

where uj£, is the driving frequency, and where s{t) is a differentiable spline in- 
terpolation function on the time interval [0,r] with vanishing first derivatives 
at both ends, which is monotonically increasing from to 1. This procedure 
results in a giant oscillation of the nucleus at t = r, provided that r is chosen 
T ^ uj]j^. The deformation is dependent on the value of the coupling kq 
chosen. 

We now use the initial conditions generated in this way (with r = 1500 
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fm/c) to study the time evolution under the action of our Hamihonian as 
defined in Eq. (|I]). The upper panel (a) of Figure 1 contains the results of 
our calculations. One can clearly observe a regular undamped oscillation 
of the quadrupole moment in coordinate space as a function of time. One 
can observe that the period of oscillation has been stretched from the 0- 
coupling value of 88.4 fm/c to 128 fm/c. This is consistent with the analytic 



calculations for infinitesimal deformations 16], which yield 



Q2{t) = Q2{to) exp{iuj2+t) 

= V2uJo (9) 

for the giant quadrupole frequency, and consequently T = 125 fm/c for 
uq = 0.0355 c/fm used in this example. 

In the lower panel (b) of Figure 1 we show our calculations for the octupole 
case. Again we see harmonic oscillations (The small variations in amplitude 
are in both cases due to beating between the initial driving frequency and 
the oscillation frequency of the self-consistent calculation.). The observed 
oscillation period is ^ 66 fm/c, in good agreement with the analytical 
result of [ig 

^^3- = V7uJo (10) 
for the giant octupole oscillation frequency, which results in T3 ^ 66.8 fm/c 



for our value of luq. 

The most important observation is here, however, that there is no damp- 
ing of the collective motion apparent in our calculation, thus indicating that 
no chaoticity is present in the collective multipole coordinates. We have 
also used slightly different initial conditions (by using a different number of 
test particles in the simulation) and obtained only slightly different results. 
This indicates that there is no sensitive dependence on the initial conditions 
present here as would be the case for chaotic motion. As an additional test, 
we performed a Fourier transform of the time signal and found one peak at 
the dominant frequency and no uj~^ noise. 

To study the single particle dynamics we employ the method of Poincare 
surface of sections. In Poincare sections, we stroboscopically record the phase 
space coordinates of particles in time. This 'stroboscope' is triggered when- 
ever the particles cross a certain plane in phase space. This technique is 
also applicable to our problem, where we solve the Vlasov equation utilizing 
the test particle method. Here we choose the plane r^; = as the trigger 
condition. 

In Figure 2 (a) we show the Poincare section for one test particle used 
to generate the octupole motion of Fig. 1 (b). One can clearly see that this 
Poincare section will become area-filling in certain regions of phase space, 
and in the limit of very long time scales. (For the present calculations we ran 
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10^ time steps of 1 fm/c each, leading to approximately 550 crossings of the 
Tx — plane from below) . In order to enable the reader to better compare the 
structures, we superimpose in Figure 2 (b) the Poincare sections of ten test 
particles with different initial conditions (but parts of the same simulation) . 
We obtain similar results for the collective quadrupole motion. 

This chaoticity is not a result of the weak destruction of the integrability 
of the corresponding static multipole potential, which can be shown to exist 
for the static octupole (but not for the static quadrupole) potential. Instead 
we attribute the chaoticity in the single particle motion to the exchange of 
energy between the motion of the individual test particles and the collective 
motion of the multipole coordinate. This exchange of energy is possible, 
because the individual test particles oscillate with frequencies, which do not 
have a rational ratio with the frequency of the collective coordinate. This 
results in the particle reaching meta-stable or unstable points in phase space 
during the course of its time evolution. At these points small changes in the 
initial conditions will have a large effect on the subsequent dynamics. An 
example for this would be the decision if the particle will temporarily oscillate 
in or out of phase with the collective coordinate. Consequently, these points 
provide large positive contributions to the Kolmogorov entropy, and chaotic 
single particle dynamics results. 

In turn one also expects each single test particle to have a randomly 
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fluctuating effect on the energy contained in the motion of the collective co- 
ordinate. However, since there are quite many test particles, these chaotic 
random fluctuations are averaged out leaving only a smooth sinusoidal os- 
cillation of the collective coordinate. This is qualitatively new in our inves- 
tigation: the generation of regular dynamics for the collective variable, the 
multipole moment of the collective oscillation, from the ensemble of single 
particles with chaotic trajectories. This is an example of how ordered macro- 
scopic motion can result form underlying chaotic microscopic dynamics. (To 
obtain this result, it was crucial to employ a self-consistent treatment of the 
dynamics entailing conservation of total energy.) 

Our finding are not restricted to the static harmonic potential, which 
we used as an example here, but they hold for a general class of central 
potentials, Vq. We have also performed calculations for a Vq oc potential 
with very similar results to the ones presented here [ED . 



And, finally, one may speculate that this interplay between chaoticity in 
individual single particle degrees of freedom and regularity in certain collec- 
tive coordinates may also play a role in the time evolution of other physical 
systems. Examples that come to mind as likely candidates are plamas in 
a tokamak, the human brain wave activity, the weather. Chaos on a mi- 
croscopic level need not necessarily lead to a catastrophic breakdown of the 
system on the macroscopic scale. 
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Figure Captions 



Figure 1 

Time evolution of quadrupole (a) and octupole (b) deformations of an 
^4=200 nucleus in coordinate space under propagation with the Hamil- 
tonian of Eq. (|l]) with a static harmonic oscillator potential and self- 
consistent coupling strength A2 = X2,s = 8.9 ■ 10~^ MeVfm""^ (a), and 
A3 = A3,, = 1.3 ■ 10-4 MeVfm-6 (b). 

Figure 2 

Poincare section of a single test particle (a) and an ensemble of 10 test 
particles (b) from the calculation leading to the time evolution of the 
octupole moment shown in Fig. 1 (b). 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9310025vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9310025vl 



